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^ \ Abstract 

Jupiter family comet 17P/Holmes underwent a remarkable outburst on UT 2007 
Oct. 24, in which the integrated brightness abruptly increased by about a factor of a 
million. We obtained near infrared (0.8 - 4.2/im) spectra of 17P/Holmes on UT 2007 
Oct. 27, 28 and 31, using the 3.0-m NASA Infrared Telescope Facility (IRTF) atop 
p '■ Mauna Kea. Two broad absorption bands were found in the reflectance spectra with 

. centers (at 2//m and 3/im, respectively) and overall shapes consistent with the presence 

of water ice grains in the coma. Synthetic mixing models of these bands suggest an 
origin in cold ice grains of micron size. Curiously, though, the expected 1.5/im band of 



Or 



water ice was not detected in our data, an observation for which we have no explanation. 
' Simultaneously, excess thermal emission in the spectra at wavelengths beyond 3.2/im 

^3 . has a color temperature of 360 ± 40 K (corresponding to a superheat factor of ~ 2.0 ± 

0.2 at 2.45 AU). This is too hot for these grains to be icy. The detection of both water 
ice spectral features and short-wavelength thermal emission suggests that the coma of 
! 17P/Holmes has two components (hot, refractory dust and cold ice grains) which are 

\ not in thermal contact. A similarity to grains ejected into the coma of 9P/Tempel 1 by 

the Deep Impact spacecraft is noted. 
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■ 1. Introduction 



X 



As a Jupiter- family comet, 17P /Holmes is unr emarkable in terms of either its dynamical prop- 



erties or its chemical composition ([Whipple! Il984h , However, on UT 2007 Oct. 24, its brightness 
increased by about a million t imes in less than a day, with its apparent magnitude rising from 



17 to ~1.5 (jBuzzi et al.l 120071 ). In comparison, NASA's Deep Impact into the nucleus of comet 



9P/Te mpel 1 in 2005 change d the brightness of the comet by a comparatively modest two magni- 



tudes ([Fernandez et al.l 120071 ). The spectacular brightening of 17P echoed an earlier outburst on 
November 6 1892, which led to the discovery of the comet by the British amateur astronomer Edwin 
Holmes in London. The two outbursts, 115 years apart, are strikingly similar in the sense that both 
occured 5 months after the comet had passed its perihelion and in both cases the cometary coma 
increased dramatically in angular size and brightness. However, whereas the comet experienced 
a second outburst in January 1893, two months after the initial event, no such second outburst 
followed the October 2007 event. 
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Over the year s, various models a nd hypotheses have been advanced to explain cometary out- 
burst phenomena (jGronkowskil 120071 ) . Some outbursts are associated with cometary splitting, in 
which the primary nucleus ejects discrete fragments, perhaps in response to rotational or other 
i nstabilities. Unfortun ately, there is no strong evidence concerning the rotational period of 17P 
(jSnodgrass et al.ll2006l ) and nothing to suggest that it is close to rotational break-up. Some comets 
split due to strong tidal forces from nearby massive objects, usually the Sun or Jupiter (IBoehnhardt 
2004 ). However, 17P/Holmes was far f rom Jupit e r and the Sun at the time of the outbursts, giving 
no reason to suspect a tidal trigger. IWhippld (| 19841 ) described a scenario whereby the impact 
of an unseen satellite could have triggered the 1892 outburst. However, the repetitive nature of 
the outbursts in this comet renders this explanation implausible. A conjectural scenario is that a 
large fraction of 17P/Holmes is amorphous water ice in which are trapped considerable amounts 
of volatile molecules. When close to the Sun, elevated internal temperatures might have precipi- 
tated an irreversible phase transition in water ice from amorphous to crystalline. This (exothermic) 
transition process would release the trapped volatiles, driving up the gas pressure and eventually 
triggering the observed outburst (e.g. Prialnik et al. 2004). 

In this paper we study the material ejected from 17P via near infrared spectroscopy (0.8 - 
4.2/im) obtained shortly after the start of the outburst. 



2. Observations 

Observations were obtained using the NASA Infrared Telescope Facility (IRTF) 3-m tele- 
scope atop Mauna Kea, H awaii. We observed 17P/Holmes on UT 2007 Oct. 27, 28 and 31, using 



SpeX (jRayner et al.ll2003l ). a medium-resolution 0.8-5.5/im spectrograph. SpeX is equipped with a 
Raytheon 1024 xl024 InSb array having a spatial scale of 0'.'15 pixel" 1 . Two cross-dispersed modes, 
known as SXD (0.8 - 2.4/im) and LXD (1.9 - 4.2/xm), were used to cover an overall wavelength 
range from 0.8^m to 4.2/im for all our IRTF observations. The spectral resolution afforded by 
SpeX is dependent, in part, on the width of the spectrograph slit. To achieve high resolution, we 
used a C/3X 15" slit that provided an average spectral resolving power of A/A A = 2300. The slit 
was projected North-South on the sky. Differential refraction is minor because 17P/Holmes is an 
extended source. To avoid potential contamination from the coma, we used a large nod distance of 
10 arcminutes to measure the sky background. 



The SpeX data were reduced using the reduction pipeline SpeXtool (jCushing et al.1 |2004j) . 
Individual flats were scaled to a common median flux level, then combined and normalized to 
generate a master flat field frame. Bad pixels were identified as outliers using flat field frames 
as well as identified from a stored bad pixel mask. Pixels with poor response were replaced by 
interpolation over neighboring normal pixels. For standard stars, the first-order sky removal was 
achieved via subtracting image pairs, 7.5 arcsec. apart, with the object dithered along the slit. 
For the comet, the first-order sky was removed by subtracting each sky (so-called B-beam) image, 
10 arcminute East of the comet, from the corresponding object (A-beam) image. The pixel-to- 



-3 - 



pixel sensitivity variations were calibrated by dividing the master flat field frame by the object 
frames. The wavelength calibration was accomplished using argon lines for the SXD mode, while 
a combination of argon lines and several sky emission lines were used for the LXD mode. Finally, 
one-dimensional spectra of the comet and standard stars were extracted from the two-dimensional 
sky-subtracted images. The SpeXtool pipeline provides two extraction modes that are designed 
for point sources and extended sources respectively. The difference between these two modes is 
that the extended-source mode (ESM) sums up all the pixel counts uniformly within the extraction 
aperture, whereas the point-source mode (PSM) specifies the aperture center first and then carries 
out a profile extraction, weighted by the spatial profile of the object along the slit. Although 
17P/Holmes was, strictly speaking, an extended source, the spatial profile of this comet is in fact 
highly centrally-condensed and appeared roughly Gaussian (Figure [1]). Therefore, we performed 
test extractions using both ESM and PSM. We found that the two modes produced similar results, 
but the PSM provided smoother spectra with slightly higher signal-to- noise ratios (SNR). Therefore, 
we adopted the point source extraction mode using an average extraction aperture box width of 
9.0", to maximize the SNR. 

In order to remove telluric absorption features we obtained spectra of a "telluri c standard star" 



close in both time and sky position to the comet. We chose the A-type star BS 1261 (jPerrvman et al 



19971 ) as the telluric-calibration standard star, which also served as a reference to approximately 
flux calibrate the comet spectra. We chose an A-type star, because such early-type stars generally 
have few metal lines and can be well-modeled by a blackbody in the near infrared (e.g. Vacca et al. 
2003). Once the target spectra had been calibrated for telluric absorptions, the different spectral 
orders from SpeX were merged to construct a single, continuous spectrum covering the wavelength 
range from 0.8 to 4.2/im. However, the spectra at the longer wavelength end suffered from the 
rapid decrease in the sensitivity of the spectrograph and appeared noisy. Therefore, we discarded 
the low quality data in the wavelength region at 4.0/mi to 4.2/iin. Since no absolute flux-calibrated 
solar spectrum is available at wavelengths longer than 2.5/im, we used a spectrum of the solar 
analog HD 76151 from the IRTF spectral library (Rayner et al., 2008)0 as a substitute for the solar 
spectrum. As a check, we carefully examined the discre pancies between the absolute flux-calibrated 



solar spectrum (0.8 - 2.5/im) from lColina et al.l (|1996l ) and the spectrum of HD 76151. We found 



that the difference between these two is less than 2% over the 0.8/im to 2.5/xm wavelength range, 
which is negligible for our present purposes. 

A journal of observations is provided in Table 1. 



http://irtfweb.ifa.hawaii.edu/~spex/WebLibrary/index.html 
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3. Analysis 

3.1. Thermal Emission 

We computed reflectance spectra of the comet via dividing the flux calibrated spectra of 17P 
by the library spectrum of the solar substitute HD 76151 (Rayner et al., 2008). The reflectance 
spectra of 17P/Holmes rise sharply beyond 3.2/iin (see Figure [2|), which we interpret as the result 
of thermal emission from hot dust in the coma. Beyond 3.2/im, the spectrum is the sum of two 
parts, one from scattered sunlight and the other from thermal emission. The scattered and thermal 
components were modeled separately from each other. Firstly, we assumed that the spectrum at ^ 
2.5/um consists of pure scattere d light and the scattering flux was estimated using an inverse-square 



law described in (|Russelllll916l ) 



fQ(X)p(X)(j)(a)C 

Jscat{* — 2A~2 " 

where f SC at{X) (W m~ 2 A~ 1 ) is the flux density of scattered solar radiation and /q(A) (W m~ 2 A~ 1 ) 
is the flux density of the solar radiation at 1AU, p(X) is the geometric albedo of the grains in the 
comet and 4>(a) is the scattering phase function. Previous studies show that the cometary phase 
function can be adequately represented by: 

</>(a) = 10~ a4/3a (2) 



wher e a (degree) is the phase angle and (5 (mag deg 1 ) is the linear phase coe fficient (IJewitt &; Meech 



19881 ). When the phase angle is small, a ^ 30°, the measured f3 ps 0.03 (jlVleech fc Jewittl Il987h . 
In Equation ([I]), r (AU) and A (m) are the heliocentric and geocentric distances, respectively, and 
C (m 2 ) is a scale factor. Physically, C is the total cross section of the cometary grains within the 
projected slit profile. Due to the uncertainty of grain sizes and albedo information, we treated C 
as a free fitting parameter. 

We estimated the slope of the reflectance spectrum of 17P using regions judged to be free 
from major atmospheric absorptions and potential water ice absorption bands, namely from 0.9 - 
1.1/xm, 1.7 - 1.8/xm and 2.2 - 2.3/iin. We made the solar analog spectrum more blue and scaled it to 
represent the scattering flux. Our scattering model, shown in Figure EJ fits the observed spectrum 
well. Assuming the slope of the comet's albedo with wavelength remains the same from 0.8fim 
to 4.2/im, we extrapolated the fitted spectral slope from the short wavelengths into the thermal 
wavelengths. Finally, we subtracted the modeled scattered light from the observed spectrum leaving 
the residual shown in Figure HI which is assumed to be pure thermal emission from the coma. Only 
the data from UT Oct. 27 were used for the thermal modeling, owing to substantial fading of the 
coma on later dates. 



W e modeled the thermal flux density, following the method described in 

Eaai): 



Jewitt &: Me cch 
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/bb(A) 



e(X)B x (T e )C e 
A 2 



(3) 



where C e is the emitting cross-section and B\(T C ) is the Planck function evaluated at the effective 
grain temperature, T e . The emissivity, e(A), depends on the ratio of the grain size to the radiation 



wavelength and o ther factors that are poorly constrained at present (jBohren fe Huf fman 



Mann et al.ll2006l ). Lacking prior knowledge of the grain sizes, we adopt e(A) = 1.0 (jGehrz et al 



1983: 



20051 ). which means that the effective grain temperature is equal to the color temperature. The 
best fitting model, shown in Figure H] by a red line, suggests that the color temperature of the 
thermal component is T c = 360 ± 40 K. As shown in Figure EJ our "scattering+thermal" model 
fits the observed comet spectrum adequately well. The considerable uncertainty on T c results from 
the fact that only a limited range of thermal wavelengths could be obtained, and that this range 
does not encompass the Planck maximum. 

For comparison, the expected isoth ermal blackbody tem perature is Tbb =180 K at r = 2.4 AU. 
The "superheat factor", S, is defined in lGehrz fc Nevl (|1992l ) to quantify the ratio between the color 
temperature of a comet and the temperature of a spherical, co-located isothermal blackbody, which 
is expressed as S = T c /Tbb- Therefore, we find S = 2.0 ± 0.2 for P/Holmes. Most comets have 



1.8 (Mason et al 



S in the range of 1 ^ S ^ 1.5 (ISitko et al,ll2004l ). but comet Hale-Bopp had S 
similar to the value measured here. 

The superheat depends on many unknown grain properties and it is not possible to use it 
to make definite statements about the grains. However, it is informative to consider the simplest 
model, in which the e missivity scales linearly w ith particle size e(A) ~ =^ as expected for particles 
in the Rayleigh limit (jBohren &: Huffmanlll983l ). Then, the temperature is related to the emissivity 
through the energy balance equation for an isothermal sphere: 



^-(l-A x )=4e x aT 4 (4) 

where F@ is the solar constant (1.36 x 10 3 W m -2 ), r (AU) is the heliocentric distance, A\ is 
the effective bond albedo at a wavelength corresponding to the peak of the solar spectrum (A ~ 
5000 A) and e\ is the effective emissivity at a wavelength corresponding to the peak of the dust 
emission spectrum (i.e. A ~ 10 /jm for T ^300 K). For a blackbody, Abb = 0, ebb = 1, we have 
Tbb = (t - ? - ) 1 ^ 4 - And so, the temperature can be written as: 



'1 - A^\ 1/4 

( i— (5) 



nvmg 
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I -A x 



1/4 



(6) 



Substituting S = 2, and with A\ -C 1, Equation © gives ~ 0.06 and so a ~ |^ ~ 0.1 //m 
as the characteristic particle size. This estimate is clearly very crude and the particle size cannot 
be trusted, probably even to within a factor of a few. However, the important point is that the 
large super heat suggests that the dominant emitters on UT Oct. 27 were small, as they were in 



Hale-Bopp (jMason et al. 



2001 



3.2. Water Ice Absorption 



Unlike the depth of an absorption feature, the band profile is relatively insensitive to continuum 
fitting and thermal emission removal. Therefore, the shape of the absorption featur es is useful as a 



constr aint on the physical characteristics of the corresponding grain s, as di s cussed in ISunshine et al 



(|2007l ) . The synthetic water ice spectrum was computed based on iHapkd ()1981 



199 



3). 



w 



Mo 



4-7T fj, + fJL 



{[l-B(g)]P(g) + H(p )H(^)-l} 



(7) 



where R is the bidirectional reflectance, fio = cos{i) and \x = cos(e), % and e are the angles of 
the incident and emergent light, respectively, g is the phase function between i and e, and w is 
the average single scattering albedo of the surface particles. Parameter B(g) is a backscattering 
function that describes the mi crostructure of the surface (e.g. porosity, grain size distribution 
and spacing b etween particles) (IDavies et al.lll997l ). P(g) is the particle phase function, H(/x) and 
H(^o) ar e the IChandrasekharl (|1960l ) H function at angles \i and /j,q, respectively. To limit the 
total number of free parameters, all the calculations are at zero phase (g = 0). In our models, we 
assumed abundant fine grains in the coma, giving B(0) ~ exp(-w; 2 /2), an isotropic phase function 
[i.e. P(0)=1] and no internal scattering within individual particles. In addition, we assume that 
the grains in the coma are spherical and smooth. Thus, the internal scattering coefficient S/ is the 
same as the external scattering coefficient S^. We used optical constants of water ice from Warren 
(1984). 

To constrain the abundance and the size distribution of icy grains, we used a "spatial mixing" 
model. The spatial mixing (also known as linear mixing) model, assumed that each component is 
physically well separated from others and so we linearly mixed the pure water ice with dark and 
featureless refractory components (e.g. amorphous carbon). As shown in Figure [6^., our water ice 
model successfully fits the two absorption features, in terms of band centers and overall shapes. In 
addition, fine structure (the Fresnel peak) was observed near 3.1/im, as seen in Figure Eb, where 
the red line is a smoothed spectrum to guide the eye. This feature is due to front-surfa ce reflection 



from ice grains (jBrown et al.l 120061 ) . and is observed in some Saturnian icy satellites (jClark et al 
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2008 ) . Unfortunately, such a detailed feature is difficult for a single-sized linear mixing model to 
match. Our best-fit model suggests that the average grain size for 17P/Holmes is about 2/zm and 
the water ice abundance is about 30% by surface area. However, one needs to bear in mind that 
the derived ice abundance is mainly determined from the depth of the 3/im band, which can be 
affected by the continuum fitting and the removal of thermal emission. Therefore, our estimate of 
the water abundance is just a lower limit to the total amount of water ice that could be in the 
coma. 

The linear mixing model adequately fits the 2 and 3/im features but it differs from the data at 
shorter wavelengths (< 1.8/im). Water ice, the major component in the synthetic model, produces 
an additional absorption feature at 1.5//m, only slightly weaker than the 2.0/im band. Furthermore, 
an additional sharp minimum at 1.65/xm, when detected, is widely used as an indicator that the 
water ice is crystalline in structure. In the spectra of I7P/Holmes, we detected a fairly strong 2/mi 
band, however, the 1.5//m band and the 1.65/im feature were absent. 

We explored two ways to decrease the strength of the 1.5/xm band. First, reducing the grain 
size can significantly diminish this feature. However, the strengths of the 1.5/tm and 2.0//m bands 
vary together and we found it impossible in our models to selectively weaken the 1.5/im band 
while preserving the one at 2//m (See Figure [7^.). Second, we explored the influence of absorbing 
impurities in suppression of the 1.5//m band. Again, these experiments affected the 1.5/um and 
2.0/um bands similarly and we were unable to match the 1.5 //m region of the spectrum of P/Holmes 
using dirty ice grains (See Figure Eb.). Conceivably, some combination of grain size and absorbing 
impurity could selectively diminish the 1.5/xm band, but the parameter space is huge and we have 
been unable to reach a solution. 

Is it possible that the 2.0/tm and 3.0/im bands in the spectrum of P/Holmes do not reflect 
the presence of water ice? We examined the spectrum of numerous other materials in order to 
find a match with the P/Holmes spectrum, without success. We closely considered the possibility 
that the bands might indicate the presence of a hydrated mineral, rather than pure water ice, since 
hydrated minerals are known to show a band near 3/iin. However, the absorption bands in hydrated 
minerals are shifted by intermolecular forces to wavelengths significantly different from the band 
centers observed in water ice and in P/Holmes. For example, the cl ay mineral m ontmorillonite 



shows a deep band starting at 2.8/iin with a minimum near 2.9/um (jRoushl 120051 ). incompatible 



with the comet data. Furthermore, the spectr a of hydrated minerals (s uch as nontronite and 



palagonite) usually show a sharp band at 1.9/Ltm (jMilliken &: Mustard! 120071 ). inconsistent with the 



2pim band in P/Holmes. These different band centers and shapes would be easily resolved in our 
data. The positions, shapes and relative depths of the 2.0/xm and 3.0//m bands point strongly to 
their origin in water ice rather than in any known hydrated mineral. This concl usion is consistent 



with in-situ meas urements of the dust in comets P/Wild 2 (jKeller et al.l 12009 ) and P/Tempel 1 



(jLisse et al.ll2007l ). where hydrated minerals are at best a trace species with a very minor impact 
on the reflected spectrum. Therefore, we proceed on the interpretation that the 2.0/tm and 3.0/mr 
bands are indicative of water ice grains, while acknowledging that we are left with the unsolved 
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mystery of why the 1.5pm band is undetected in P/Holmes. 

Although the comet was faint for the LXD observation, we were able to obtain SXD data of 
high quality on UT Oct. 28 and 31, 2007. We observed the 2.0pm water ice absorption feature in 
all the SXD data, shown in Figure EB The relatively long survival time of this ice feature at the 
heliocentric distance of 2.4AU provides an extra constraint on the physical properties of the coma 
grains. To estimate the lifetime of cometary grains, we calculated the sublimation rate as a function 
of grain size by assuming an equilibrium state between heating (due to absorption of solar radiation) 
and cooling (due to re-radiation and sublimation) at the surface of the grain: Eq = E ra d + E su b. 
Here, Eq (W m -2 ) is the heating rate per unit area, given by: 



E e = l°° ^pr-Q^ dX ( 8 ) 

where /©(A) (W m~ 2 A^ 1 ) describes the solar flux density, r(AU) is the heliocentric distance and 
Qabs is the absorption efficiency, which was computed using Mie theory. The quantities E ra d (W 
m~ 2 ) and E su j, (W m~ 2 ) are the re-radiation rate per unit area and the sublimation rate per unit 
area respectively, which can be expressed as: 



E 



rad 



B(X,T s )Q e dX 



(9) 



and 



E 



sub 



HP(T S 



m 
2irkT s 



(10) 



where H (J kg 1 ) is the latent heat of sublimation, P(T S ) (N m 2 ) is the equilibrium vapor pressure 



(jFanale fc Salvaillll984l ) at the surface temperature T s , m (kg) is the mass of a water molecule and 
k (J K _1 ) is the Boltzmann constant. Also, B(X,T S ) is the blackbody radiation at the surface tem- 
perature T s and Q e the emission efficiency. Under the assumption of thermo dynamic equilibrium 
betw een the particle and the surrounding radiation field, we have Q a bs = Qe (jBohren fc Huffman 
19831 ). Then, equations ([8))- (|10p can be compiled in the form: 



HP(T S 



m 



27T/CT, 



/o(A) 
4r 2 



B(X,T S 



Qabs dX 



(11) 



We solve this equation to find the grain temperature T s and then compute the specific sublimation 
rate ip (kg m~ 2 s" 1 ), where ip = P(T S )^J 27 ^t • Consequently, we can calculate the lifetime of an 

ice grain by r = p J^ M r 4)~ 1 da, where p is the water ice density and a m , au are the minimum and 
maximum grain diameter, respectively. 



- 9 - 



For pure ice grains, the sublimation rate was calculated using optical constants from Warren 
(1984). However, the observed cometary grains may not be pure but contaminated by refractory 
particles, such as amorphous carbon. Impurities can significantly shorten the lifetime of an icy 
grain by reducing the albedo and in creasing the equilibrium t emperature in sunlight. The optical 



constants of amorphous carbon from lRouleau fc Marti n (1991 ) were used as the darkening material 



in the dirty ice and the Maxwell- Gar nett mixing rule (|Mukailll986l ) was applied. The sublimation 



lifetimes of pure and dirty water ice grains as functions of their size are shown in Figure El 

First of all, the presence and strength of the water ice absorption features indicates that the 
fresh icy grains survived against sublimation over a substantial fraction of the projected area of 
the slit. We extracted spectra of pure coma from a region that is l'/5 to 3" away from the nucleus. 
Although the SNR of the outer coma became much weaker, we found the water ice features persisted 
in the coma spectra. Our result suggests that the ice grains survived the sublimation up to 3" away 
from the nucleus. The geocentric distance was 1.62 AU at the time of the observation, so that the 
projected distance of 3" on the coma corresponds to 3534 km. The measured expansion speed of 
the coma is about 500 m/sec (Rachel Stevenson 2009, private communication), and so it would 
take about 2 hours for newly released icy grains to travel from the nucleus to the region that is 3.5 
thousand kilometers away, if moving parallel to the length of the slit. As noted above, the best-fit 
to the band shapes at 2/im and 3/Lim gives an average diameter of coma ice particles of ~ 2^m. 
Figure [8] shows that dirty ice grains of this size are not able to survive for 2 hours, the residence 
time in the slit. More importantly, the source of the observed ice grains likely survived for at least 
7 days after the initial outburst. The observed long-lasting water ice absorption feature suggests 
that the ice grains in the coma of 17P must be relatively clean. 



3.3. Spectral Slope 

The negative slope (Figure [9]) distinguishes 17P from most comets that gene rally show neutral 



or positive slopes at optical and near infrared wavelengths (jJewitt k, Meechll 19861 ). Previous studies 
have shown that near infrared spectra of pure water ice or hydrated minerals also have negatively 
sloped continua within the same wavelength range (Clark et al. 1986). The spectral slopes in 
P /Holmes grew steeper with time from the initial outburst, as shown in Figure suggesting that 
either the abundance of icy grains was increasing or the average grain size was decreasing between 
the nights. Since the observations were made after the outburst and the gas outflows became 
weaker with time, the shrinkage of grains is more likely to be responsible for the change of the 
spectral slope. The scattering efficiency, Q s , generally depends on the ratio of the particle size to 
the wavelength. When the grain sizes are comparable to or slightly smaller than the wavelength, the 
scattering efficiency varies inversely with the wavelength. Therefore, as the average size of coma 
particles becomes smaller, the blue scattering efficiency increases and the reflectance spectrum 
becomes steeper. 
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3.4. Other Spectral Features 



In addition to water ice absorptions, a broad absorption feature was found near 1.2/mi in the 
spectrum of 17P/Holmes. To verify the presence of this absorption feature, we examined the SXD 
data of 17P/Holmes obtained on UT 2007, Oct. 27, Oct. 28 and Oct. 31, using two different solar 
analog stars and the flux calibrated solar spectrum to calculate the relative reflectance. We found 
that the 1.2//m feature appeared consistently over three nights and its existence was independent 
to the choice of the solar analog. Furthermore, the center and the overall profile of this absorption 
feature are also insensitive to the particular solar analogs used. Some alkanes, such as methane 
and ethane, show an absorption band near 1.2/im attributed to the second overtone of C-H stretch. 
However, these hydrocarbons also show other, stronger features near 1.7 and 2.3/xm that do not ap- 
pear in the spectrum of 17P/Holmes (see Figure [T0|) . So, we do not favor an interpretation in terms 
of hydr ocarbons. Fu rthermore, we checked for this feature amongst the diffuse interstellar bands 
(DIBs) (|Herbiglll995l ). but found no matches. At present, the 1.2//m feature remains unidentified. 



Significant emission features were detected near 3.4/um in spectra of 17P/Holmes from UT 
2007 Oct. 27 and 28. Although the data of UT 2007 Oct. 28 appeared much more noisy than from 
UT 2007 Oct. 27, the 3.4/im emission feature was strong enough to be discerned clearly from the 
noisy continuum. The LXD data of UT 2007 Oct. 31 suffered from very low SNR and thus were not 
used. As shown in Figure [TTT the peak wavelength of the stronger emission band is near ~ 3.36/mi 
and the width is ~ 0.15/zm. In addition, there is a weaker feature centered at 3.52^m with the 
width of 0.07/zm. These characteristics are consistent with the ^2, ^3 and z^g CH-stretching bands 
of methanol , which were first observed i n comet Halley two decade s ago and seen in several comets 

). W e calculated the methanol 



since then (Hoban et al 



1991 



Reuter 



1992; 



Mumma et al 



2001 



production rate following the method described in iDisanti et al.1 (|1995l ). The derived production 
rate is about 2.2 ± 0.5 (10 27 mol/s), which is comparab le with those of o t her co mets. However, 
we note that the methanol production rate measured in Dello Russo et al. (j2008l ) is several times 
higher than the value we obtained in this work. This difference may due to the timing difference 
between their observations and ours. Moreover, the measurements of Dello Russo et al. are based 
on high-resolution observations. For these reasons, we do not see a significant contradiction between 
their results and ours. 



4. Discussion 

Our thermal model suggests that the color temperature of the coma of P/Holmes is about 
360K, however, water ice absorption features were also detected in the spectra at the same time. 
Given that the life time of ice grains decreases as the temperature increases, the presence of ice 
grains requires a low temperature (T < 270 K) to ensure the survival of solid water. One explana- 
tion for this temperature dichotomy is that the coma of P/Holmes might contain two distinctive 
grain populations that have different temperatures not in thermal equilibrium at the time of the 
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observations. The temperature of a particle at a given heliocentric distance is determined by many 
factors, including the chemical composition and the grain size. For sufficiently small absorbing 
grains (e.g. c arbonaceous compounds), the size exerts a dominant influence on the temperature 
(|Hannerlll983l ). Most particles made of transparent mate rials have equilibrium t emperatures below 
that of a co- located blackbody regardless of their size (jKolokolova et al.l 12004 ) . Compared with 
carbon, pure water ice has much higher heat capacity and albedo. If water ice grains and refrac- 
tory grains are not in thermal contact, they could maintain distinct temperatures. Therefore, our 
observations indicate that the coma of 17P/Holmes consists of distinct cold, icy grains and hot fine 
refractory grains. 

Interestingly, our inference of the coexistence of cold and hot particles in the coma of 17P/Holmes 
is consistent with observations of P/Tempel 1 from NASA's "Deep Impact" (DI) mission. There, 
a strong water ice absorption at 3/xm was observed in the spectra of the ejecta within 3 sec of 
the i mpact. Synthetic mix ing models revealed small, pure particles with an average diameter of 
lfim ([Sunshine et al.l 120071 ). Simultaneously, a significant thermal component from hot dust was 
observed along with the ice features in the DI spectra at wavelengths > 3/xm. Therefore, the impact 
liberated ejecta is likely to contain both pure water ice grains and hot dust particles. In addition, 
the survival of water ice features throughout the DI flyby observations indicates th at the ice grains 



were thermally /physically separated from other components and therefore pure ([Sunshine et al 
20071 ). 



In other comets, fragmentation of grains is believed to be an important process. For exam- 
ple, the distri bution of impacts of small dust grains onto s pacecraft in the comae of P/Halley 



and P/Wild 2 (jBoehnhardt & Fechtie 



1987 



Clark et al.l 12004 ) was highly non-random, suggesting 



spatially correlated swarms of particles that could have been produced by fragmentation of larger 
parent particles. Fragmentation could result from electrostatic disruption caused by the build-up 
of charge on aggregate dust grains exposed to the ionising radiation of the Sun. Alternatively, 
fragmentation could be due to the loss of a binding agent ( "glue" ) holding sub-grains together in 
the initial aggregate. 

These considerations suggest a simple picture in which the body of 17P/Holmes consists of 
porous aggregates of refractory particles and ice grains held together by ice and / or organic com- 
pound s (the latter suggested by the detection of the 3.4//m methanol feature; see also lLamv fc Perrin 
(jl988l )). Once ejected from the nucleus into sunlight, the ice begins to sublimate and the aggregates 
break up, releasing the particles we detected and contributing to the dramatic increase in brightness 
of the coma. A porous aggregate structure of low tensile strength could be representative of the 
nucleus as a whole, and fragmentation of ejected pieces of the nucleus might explain the multiple 
secondary sources observed in high-contrast imaging of the coma. 
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5. Summary 

Near infrared spectral observations of comet 17P /Holmes were obtained within a week of its 
UT 2007 Oct. 24 outburst, with the following results: 

• Two deep absorption features were observed at 2^m and 3(im respectively, consistent in 
central wavelength, band width and shape with the presence in the coma of water ice grains. 
Spectral fits with linear mixing models show that 17P/Holmes contains a significant fraction 
of pure water ice and that the central coma is dominated by fine ice grains of micron size. 

• The spectral fits with water ice are so good that we find it difficult to believe that the 2/im 
and 3^m bands could have any other explanation. Indeed, an explicit search for hydrated 
and other minerals that have similar features failed to produce a convincing match. However, 
with water ice models, we have yet to explain why the 1.5//m band of water ice is absent from 
our data. 

• We observe strong thermal emission at wavelengths A > 3.2//m, which we attribute to hot 
dust. The derived color temperature is 360 =L 40K and the corresponding superheat factor is 
~ 2.0. As in other comets (e.g. C/Hale-Bopp) the high superheat indicates that the thermal 
emission is dominated by small, refractory grains. 

• The simultaneous presence of hot, refractory grains and water ice grains (which must be 
cold in order to minimize sublimation) shows that these two populations are not in thermal 
contact. 

• A solid-state absorption band observed near 1.2/im on three nights remains unidentified. 
Further laboratory experiments are needed. 

• Our observations of 17P/Holmes resemble those of P/Tempel 1 from NASA's Deep Impact 
mission, in the sense that both comets showed small, pure ice grains coexisting with hot 
refractory grains in the coma. 
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Table 1. Observational parameters 



UT Date 


Wavelength 


Itime 


Airmass 


r 


A 


a 


CSO Tau 




fim 


(s) 




AU 


Ai: 


deg 




2007 Oct. 27 


0.8-4.2 


1440 (SXD) + 780 (LXD) 


1.16 - 1.32 


2.44 


1.63 


16.2 


0.05 


2007 Oct. 28 


0.8-4.2 


1200 (SXD) + 900 (LXD) 


1.29 - 2.30 


2.45 


1.62 


15.9 


0.04 


2007 Oct. 31 


0.8-4.2 


720 (SXD) + 450 (LXD) 


1.16 - 1.19 


2.46 


1.62 


15.1 


0.06 
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t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




j i I i i I I i i i I i i i I i i i I i i i I i i i I i L 



-6 -4 -2 2 4 6 
Spatial Profile (arcsec) 

Fig. 1. — The solid line shows the spatial profile of the comet along the slit. The two dashed lines 
show the regio, of width of 9.0", from which the final spectrum was extracted. Although the comet 
was surround by a sizable coma at the time of the observation, the surface brightness profile of 17P 
was centrally condensed, with a roughly Gaussian profile having a FWHM of 1.6". 
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Fig. 2. — The reflectance spectrum of 17P/Holmes from 0.8- to 4.0-/xm, taken on UT 2007 Oct. 27 
and normalized at A = 2.2/im. At shorter wavelengths (A < 3.0/um), the most significant spectral 
feature of Holmes is the negatively sloped continuum. A linear fit to the continuum is shown as 
the red dashed line. The rising spectrum at wavelengths (A > 3.2/xm) is attributed to thermal 
emission. The two green shaded rectangles block out the regions that are heavily contaminated by 
the telluric absorptions. 
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1.0 1.5 2.0 2.5 

Wavelength (micron) 



Fig. 3. — The black line shows the comet spectrum, taken on UT 2007 Oct. 27. The red line 
represents the modeled scattering spectrum, which matches well with the observed data in the 
wavelength range 1.0 - 2.5/um. 
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3.6 3.7 3.8 3.9 4.0 

Wavelength (micron) 

Fig. 4. — The black line represents the thermal component that obtained after subtraction of 
a scattering model of the spectrum of 17P/Holmes, which was taken on UT 2007 Oct. 27. The 
colored lines are spectra of blackbodies with different temperatures. The wavelength range (2.8 
< A < 3.55/im) contains a strong water ice absorption band and some emission features, so we 
only modeled wavelengths > 3.55 /jm. The best fit model, shown in red, indicates that the color 
temperature of P/Holmes is 360K. 



-21 - 




2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 

Wavelength (micron) 

Fig. 5. — The solid line shows the comet spectrum, taken on UT 2007 Oct. 27. We do not show 
the comet spectrum in the wavelength region 2.6//m - 2.8/um since it was heavily contaminated by 
the telluric absorptions. The dash-dot line is the modeled scattering flux. The long dashed line 
shows the scattering+thermal model that fits the data adequately well. Here the thermal model 
was calculated using the temperature of 360K. 
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Fig. 6. — a). UPPER A synthetic water-ice spectrum from the linear-mixing model with 2/xm 
grains is able to fit the data successfully. Spectra computed using bigger grains tend to have wider 
and deeper absorption bands than the ones were observed. The arrow points out a small peak near 
3.1/im, which is due to Fresnel reflection at the surface of ice grains, b). LOWER The black open 
squares represent the obtained spectrum of the comet on UT 2007 Oct. 27. The red dashed line is 
a smoothed curve with a smoothing width of 0.1/im, which clearly shows the Fresnel peak that is 
about 0.15/xm wide and centered near 3.05/um. 
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Fig. 7. — The black line shows the spectrum of the comet on UT 2007 Oct. 27 and the colored 
lines show linear mixing models, a) UPPER Effect of grain size, b) LOWER Effect of impurities. 
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Fig. 8. — Lifetimes of pure and dirty ice grains at the heliocentric distance of 2.45 AU. The two 
"star" symbols show the lifetimes of dirty (lower) and clean (upper) ice grains with a diameter of 
2.0/um. The lower horizontal dashed line indicates the required in-slit residence time which sets the 
lower limit of the lifetime for ice grains in the coma of 17P. It shows that only pure water ice was 
able to survive and to produce the detected absorption features. The upper horizontal dashed line 
represents the time interval between the start of the outburst and the end of our observations. 




Fig. 9. — Red, green and blue points represent the SXD spectra of 17P/Holmes taken on UT 
2007, Oct. 27, 28 and 31, which are normalized to unity at 2.2//m. Dashed lines are linear spectral 
slopes, which reflect either the chemical composition or the size distribution of the cometary coma. 
The figure illustrates the spectral slope becomes steeper with time further away from the original 
outburst. 
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Co mpar ison of the NIR spectrum of 17P/Holmes, Pluto and Charon (from 



Brownfe Calvinl (12000 )). The spectrum of Pluto is dominated by methane and water ice domi- 



nates the spectrum of Charon. The significant absorption feature near 1.2/xm in the spectrum of 
the comet is not associated with either water ice or methane. 
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Fig. 11. — Points are observed spectra of 17P/Holmes and solid curves are best fits to the contin- 
uum. On both nights, the emission features were distinct from the continuum, with a strong peak 
near 3.36/im and a weak peak at 3.52^m. 



